Synchnonization, zero-resistance states and rotating Wigner crystal 
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We show that rotational angles of electrons moving in two dimensions (2D) in a perpendicular 
magnetic field can be synchronized by an external microwave field which frequency is close to the 
Larmor frequency. The synchronization eliminates collisions between electrons and thus creates 
a regime with zero diffusion corresponding to the zero-resistance states observed in experiments 
with high mobility 2D electron gas (2DEG). For long range Coulomb interactions electrons form a 
rotating hexagonal Wigner crystal. Possible relevance of this effect for planetary rings is discussed. 
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The discovery of microwave-induced resistance oscil- 
lations (MIRO) Q and of striking zero-resistance states 
(ZRS) of a 2DEG in a magnetic field @, |j| attracted a 
great interest of the community. A variety of theoreti- 
cal explanations has been pushed forward to explain the 
appearance of ZRS (see Refs. in Q). These approaches 
provide certain MIRO which at large microwave power 
even produce a current inversion. However, these theo- 
ries do not give zero resistance, and it is usually argued 
that ZRS are created as a result of some additional insta- 
bilities which mysteriously compensate currents to zero. 
Hence, a physical origin of ZRS still remains a puzzling 
problem. 

In this work we suggest a generic physical mecha- 
nism which leads to a suppression of electron-electron 
collisions and creates ZRS. Its main element is the syn- 
chronization phenomenon which has abundant manifes- 
tations in science, nature, engineering and social life 
A simple picture of the effect is the following: a mi- 
crowave field excites electrons above the Fermi level and 
switches on dissipation processes in energy, which com- 
pensate microwave-induced energy growth, thus creating 
a noncquilibrium steady-state distribution. Due to this 
dissipation, when the microwave frequency ui is close to 
a resonance with the Larmor frequency wg, the synchro- 
nization of the phases of Larmor rotations of electrons 
with the phase of microwave field is established. In this 
way all electrons start to oscillate in phase - like male 
fireflies blink in phase on trees in Siam, emitting rhyth- 
mic light pulses in order to attract females [B|, Q ■ But 
compared to fireflies stationary sitting on trees, the syn- 
chronization of moving electrons brings a new element 
not presented in the common synchronization science: 
due to synchrony the collisions between electrons extinct 
what leads to a drastic drop of the collision-induced dif- 
fusion constant D and to creation of ZRS (we note that 
D is proportional to experimentally measured resistance 
R xx since R xy » R xx 0,0]). A simple image of such 
synchronized electrons is given by an ensemble of parti- 
cles randomly distributed on a 2D plane, which rotates 



as a whole on a Larmor circle of radius = vf/lob 
with frequency ojb- Indeed, in such a rotating ensemble 
(plane) particles never collide, and we demonstrate be- 
low that this can happen with 2D electrons synchronized 
with a microwave field phase in a magnetic field B. The 
synchronization origin of ZRS allows one also to under- 
stand qualitatively why ZRS exist only in high mobility 
samples. Indeed, it is well known that synchronization 
remains robust to a weak noise but disappears at strong 
one [H j hence a weak impurity scattering will not destroy 
ZRS. It is also important to note that the above picture 
is based on classical dynamics that has its grounds since 
in the experiments 0, Q the Landau quantum level is 
large nr, ~ 100. Thus, we start our analysis with a clas- 
sical mechanics treatment and will turn to a discussion 
of quantum effects later. 

To justify the synchronization picture of ZRS described 
above we perform extensive numerical simulations using 
two main models of classical electrons (particles) with 
short range and Coulomb interactions. In the simplest 
setup, we model particle dynamics with short range inter- 
actions in magnetic and microwave fields with the Nose- 
Hoover (NH) thermostat (see e.g. 0, HI) combined with 
interactions treated in the frame of the mesoscopic multi- 
particle collision model (MMPCM) 0. The NH ther- 
mostat produces an effective friction 7 which keeps the 
average kinetic energy (p 2 /2m) equal to a given ther- 
mostat temperature T and equilibrates heating induced 
by a microwave field f ac = f cos ujt. At the same time 
the MMPCM drives system toward ergodic state with 
the equilibrium Maxwell distribution at a given temper- 
ature. In this way the particle dynamics is described by 
the equations: 



Pi/m , p l = Fj + i Li + f ac 
7 =[(p 2 )/(2mT)-l]/r 2 
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where , Pi are the coordinate and the momentum of 
particle i, fn = e[p i x'B]/mc is the Lorentz force, Fj is 
an effective force produced by particles collisions, r is the 
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relaxation time in the NH thermostat and (p 2 ) means av- 
erage over all N particles. We usually consider the case 
of a linearly polarized microwave field f ac since numerical 
data give no significant dependence on polarization. In 
numerical simulations N particles are placed randomly 
on a square cell L x L which is periodically continued 
all over the plane. The collisions are treated in the MM- 
PCM formalism, namely the main cell is divided into N c 
small collision cells in which after a time step At the ve- 
locities of particles are reshuffled randomly but keeping 
conserved the momentum and energy of particles in the 
collision cell Q. In absence of microwave radiation the 
system evolves to a usual thermal equilibrium with the 
Maxwell distribution. The average rate D of particles 
diffusion in space is computed via their displacements af- 
ter a large time interval t. In presence of the microwave 
field the diffusion rate D is drastically changed in the 
vicinity of the resonance ujb ~ w as it is shown in Fig. [1] 
for typical values of parameters. 

Fig. [T] clearly shows the existence of a synchroniza- 
tion Arnold tongue inside which the diffusion drops to 
zero (here as well as in numerical simulations below, its 
residual numerical value D/Dq < 10~ 8 is essentially de- 
termined by roundoff errors and fit accuracy and is non 
distinguishable from zero). According to Fig. [T] the syn- 
chronization regime and the ZRS exist inside the detun- 
ing range 

\uj b -u>\< sf/(mv T ) , (3) 

where vt = y/2T/m and a numerical constant s w 0.7. 
We note that s is not sensitive to the relaxation time 
r which has been varied by an order of magnitude. In 
fact the domain of ZRS given by (O is very similar to a 
usual synchronization domain for one particle [f| which 
is also not sensitive to the dissipation rate. The ori- 
gin of this similarity is rather clear: the synchroniza- 
tion with the microwave field phase eliminates collisions 
between particles, so that they move independently and 
hence the Arnold tongue becomes the same as for one- 
particle synchronization. The fact that in the ZRS the 
collisions are eliminated, is confirmed by direct count- 
ing of the number of collisions in the numerical code 
and by computation of the synchronization parameter 

s = Ei<j( v i ~ v o) 2 /( n2v t/ 2 ) which in the ZRS dr °P s 
down to S ~ 10~ 10 being determined by roundoff errors. 
This means that all particles have the Larmor phase syn- 
chronized with the microwave field phase while their po- 
sitions in the coordinate space are disordered. Outside 
of the ZRS particles continue to diffuse with a rate D 
which is comparable with the unperturbed rate Dq. At 
small values of N c and At when the collision rate be- 
comes rather large and Dq ~ D c = v T /u> « vprs, the 
ZRS regime is destroyed. Another model of collisions, 
in which the velocities of colliding particles are changed 
randomly in a bounded relatively small scattering angle, 



gives essentially the same result ([3]) for the ZRS. 
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FIG. 1: (color online) Gray-scale plot of the normalized diffu- 
sion rate D/Dq as a function of frequency detuning (ujb— oj)/u 
and rescaled microwave field strength vjul with v — f/mvr 
where vt = \/2T/m is the thermal velocity and Do is the dif- 
fusion rate in absence of microwave at = w. The system 
parameters are: N = 1000 iV c = 4 x 10 4 , coAt = 0.2, ujt = 10, 
ut = 500, L/r B = 10, D /D c = 0.12 (with D c = v^/lu, 
p — N/L 2 and tb taken at ujb = w, thus a number of parti- 
cles inside a Larmor circle is Nb = i^r%p = Trpv^/io 2 — 10tt, 
io = const). Color intensity is proportional to D/Dq (black for 
maximum D/Do ~ 1.2 and white for minimum D/Dq — 0). 

To check the existence of the ZRS in the case of long 
range Coulomb interactions we use the molecular dynam- 
ics (MD) simulations of a classical two-dimensional elec- 
tron liquid as described in [n|. The results obtained in 
[l(| show that such an approach correctly describes plas- 
mon modes in presence of a magnetic field even when 
the Coulomb energy Ec = e 2 /a is large compared to 
classical temperature T. Here, a = 1/y/Wp is an average 
distance between electrons determined by the electron 
density p. We ensured that our numerical code with the 
Ewald resummation technique reproduces correctly the 
results presented in [10(. To equilibrate the heating in- 
duced by the microwave field we introduce in Eq. {TJ an 
energy-dependent dissipation with 7 = "/o(E — Ep)/ Ep 
for E = p 2 /2m > E F and 7 = for E < E F . In such 
a way the dynamics remains Hamiltonian for E < Ep 
while above Ep the dissipative processes are switched on 
as it is usually the case for 2DEG; thus Ep plays a role 
of Fermi energy [ll|. Usually we use Ep /T w 2 but the 
obtained results are not sensitive to this ratio. The main 
part of simulations is done at an intermediate interaction 
strength r s — Ec/Ep = 0.3 but we ensured that an in- 
crease(decreasc) of r s by a factor 7(3) does not change 
qualitatively the results (samples studied in 0, 0] have 
r s ~ 2). Also a variation of the dissipation rate 70 by 
an order of magnitude does not affect significantly the 
results and we present data at 70 ~ 0.7vp/a. The same 
is true for the total number of electrons varied from 20 
to 200 at p = const, thus we present data at N = 100. 

A typical example of the dependence of average elec- 
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FIG. 2: (color online) Dependence of electron square displace- 
ment R 2 , rescaled by electron density p, on the rescaled time 
uit. Here the Larmor frequency is ujb — oj at microwave field 
strength / = (red top curve); f /(mvFUi) = 0.059 (fa/EF = 
0.02) for ojb = w (blue bottom curve), ujb = 0.875lj (second 
from top black dashed curve), and ujb = (J with impurity 
scattering mean free path U — 96vb (second from bottom 
green curve). Total number of electrons is N = 100 and 
Nb = Tipvp/uj 2 = 34.7. The linear fit gives the diffusion rates 
D/D c = 0.089,0.068,0.0040,9 x 10~ 6 with D c = v%/u (re- 
spectively for curves from top to bottom ordered at cut — 400). 



tron square displacement R 2 on time is shown in Fig. [2] 
The introduction of microwave field leads to the synchro- 
nization of electron Larmor phases and to a drastic drop 
of diffusion rate at ujb = oj, formally by 4 orders of mag- 
nitude; the synchronization parameter S drops down to 
S « 10~ n in this case that means that collisions are 
completely switched off. A shift in the Larmor frequency 
lob = 0.875w destroys synchronization and diffusion D 
is restored being close to its unperturbed value Dq at 
/ = 0. An introduction of a weak noise linked to impu- 
rity scattering with a scattering time ti and mean free 
path li — vpTi leads to a finite diffusion rate D which is 
however much smaller than Dq until 1+ 3> tb (sec Fig.[2|. 
A decrease of the mean free path down to U « 10re de- 
stroys synchronization and restores a diffusion with rate 
D as Dq. We note that li as lOOre approximately corre- 
sponds to experimental conditions in 0, 0| . 

The dependence of D on the frequency detuning is 
shown in Fig. [3] The numerical data for Coulomb inter- 
actions between electrons show that the ZRS exist inside 
the synchronization window near the resonance lob ~ w 
with the width given by Eq. ([3]) where vt should be re- 
placed by vf and s as 0.8. Inside the ZRS the diffusion 
drops practically to zero as discussed above. The valid- 
ity of the relation ([3]) shows that the effect is not very 
sensitive to the type of interactions between particles. 

However, the long range nature of Coulomb interac- 
tions significantly modifies the structure of the ZRS con- 
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FIG. 3: (color online) Dependence of rescaled diffusion rate 
D/Do on the rescaled frequency difference (ujb — uj)/v. Here 
v — f/mvF, Do is diffusion rate in absence of microwave at 
ujb = u>, fa/EF = 0.02 and number of electrons in a Larmor 
circle is Nb ~ 2 (stars), 8 (triangles), 34.7 (squares), 138.8 
(points) with D /D c = 0.054, 0.089, 0.12, 0.14 and D /v F a = 
0.20, 0.35, 0.53, 0.64 respectively. Total number of electrons 
is TV" = 100, L = y/N/p w 17.72a. 



figuration: for short range interactions particles are dis- 
tributed over the plane in a disordered way, while for the 
Coulomb interactions electrons form a hexagonal Wigner 
crystal as it is shown in Fig. [5] The whole crystal (as 
well as each electron) is rotating in the plane with the 
frequency u> w lob and rotation radius tb = vf/ujb- 
A remarkable property of the rotating Wigner crystal is 
that formally it is formed at a rather small parameter 
r s ~ 0.3 while the usual Wigner crystal requires r s val- 
ues by more than two orders of magnitude larger 
We attribute this to synchronization of electron Larmor 
phases with the microwave field phase, what eliminates 
collisions between electrons and suppresses fluctuations, 
thus yielding an effectively large r s in the rotating frame. 
In the crystal all Coulomb forces acting on an electron are 
compensated, thus the size of synchronization domain in 
frequency range given by Eq. ([3]) is essentially the same 
as for one-particle synchronization and is practically in- 
dependent of dissipation rate 70 [l| . 

In conclusion, we have suggested a generic mechanism 
which for 2D particle rotational dynamics (e.g the Lar- 
mor rotation, but it may be also relevant for other 2D sys- 
tems like planetary rings) produces synchronization of ro- 
tational angles of all particles with the phase of external 
driving periodic field. As a result a rotating Wigner crys- 
tal is created and a collisional diffusion is suppressed by 
several orders of magnitude. The collective crystal struc- 
ture also suppresses the diffusion due to impurities. We 
propose that this effect explains the appearance of ZRS in 
2DEG observed in [1,0]. According to Eq. (0 the relative 
size of ZRS plateau is Auj/oj m 2v jui m fvp/ujEp that for 
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FIG. 4: (color online) Instant image of the rotating Wigner 
crystal formed by N = 100 electrons (points) in a periodic cell 
with L = a/A7p « 17.72a, ut = 480, lob = u, fa/E F = 0.02 
and Nb = 34.7 (as in Fig. [2] bottom curve); the circle shows 
an orbit of one electron for 240 < cut < 480; lines are drawn 
to adapt an eye showing a hexagonal crystal with a defect. 



experiments with E F ~ \QQK°, v F ~ 3 x 10 7 cm/s 

and oj/2it = ZhGHz gives Alo/lo w 0.1 if the field 
strength acting on an electron is f/e ~ 5V/ cm. This 
relative width is in a reasonable agreement with the ex- 
perimental results 0, H, EH where unfortunately an ex- 
act value of / is not known. The synchronization energy 
scale Eg ~ frg ~ lOiv and the crystal Coulomb energy 
Eq ~ 200if might be the origin of large energy scale 
EU ~ 10if in the ZRS activated transport @, An im- 
portant discrepancy from the experiments is that our the- 
ory gives synchronization only near the main resonance 
wb/u w 1 while in the experiments ZRS exist also near 
integer low resonances j with lob/oj ~ 1/j. We suppose 
that these resonances may appear due to an additional ef- 
fective 2DEG potential and surface modulation in space. 
This may generate higher harmonics of the Larmor mo- 
tion and produce synchronization also at integers j > 1. 
Such a modulation apparently appears during molecu- 
lar epitaxial growth [1J] . This modulation also produces 
a frequency shift in the rotational frequency that may 
be responsible for a resonance shift of the ZRS domain 
compared to the Larmor resonance (see H ! [Dl). A 
coherent rotation of electrons in the crystal generates a 
rotating magnetic field By/ ~ [i^evpp ~ 1G parallel to 
2DEG which can be detected experimentally. 

Our theory is based on the classical dynamics and it 
is crucial to analyze the contribution of quantum effects. 
In principle it is known that at small effective values of 
Planck's constant H e ff the synchronization is preserved 
while at large values H e f / it is destroyed by quantum fluc- 
tuations [15j. For 2DEG ft e // ~ 1/nz, and at nj, ~ 100 
it is natural to expect that the synchronization is robust 
against quantum fluctuations. However, a reduction of 
til by an order of magnitude due to an increase of cu 
to a THz range may significantly enhance quantum noise 



and destroy ZRS. Further theory development is required 
to study quantum effects properly. The most important 
question is about the amount of electrons which are in- 
volved in the rotating Wigner crystal. Indeed, our clas- 
sical studies show that all electrons are involved in this 
state but in the quantum case it is rather possible that 
only a finite fraction of electrons near the Fermi level con- 
tributes to the rotating crystal, while all other electrons 
will stay as a non-interacting background. 

Finally we make a conjecture that the mechanism de- 
scribed here may be responsible for enormously long life 
time ( ~ 10 12 rotations) and sharp edges of planetary 
rings (e.g. ~ 10m for Saturn) 1|| 17|. Indeed, a temper- 
ature there is very low and in the rotational frame the 
2D dynamics of particles is similar to motion of electrons 
in a magnetic field (lij . Hence, moons inside a ring and 
near to a resonance may produce synchronization and 
diffusion suppression with emergence of ZRS in space. 
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